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ABSTRACT 


In order to enhance the performance of pulse detonation combustors (PDCs), an 
efficient deflagration-to-detonation transition (DDT) process is critical to maintain the 
thermodynamic benefits of detonation-based combustion systems and enable their use as 
future propulsion or power generation systems. The DDT process results in the 
generation of detonation and can occur independently, but the required length is 
excessive in many applications and also limits the frequency of repeatability. 
Historically, obstacles have been used to reduce the required distance for DDT, but often 
result in a significant total pressure loss that lessens the delivered efficiency advantages 
of PDCs. This thesis evaluated various swept-ramp obstacle configurations to accelerate 
DDT in a single event PDC. Computer simulations were used to investigate the three- 
dimensional disturbances caused by various swept-ramp configurations. Experimental 
tests were conducted using various configurations that measured combustion shockwave 
speed and flame front interactions with the swept-ramp obstacles. Detonation was 
verified across the instrumented section through high-frequency pressure transducers, and 
experimental data proved that swept-ramp obstacles successfully accelerate the DDT 
process with minimal pressure losses. 
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I. 


INTRODUCTION 


Current propulsion technologies such as turbojet/turbofans and rockets have 
reached a maturity level where future development will not likely produce significant 
gains in thermodynamic efficiency and performance. Various propulsion technologies are 
capable of high performance and practical thermodynamic efficiency in a specific range 
of flight Mach number. None of the current air breathing propulsion systems have the 
capability to maintain their high performance level and efficiency over a large range of 
speed. Pulse Detonation Combustion is a technology being explored in advanced future 
propulsion concepts that will allow higher performance due to the increase in the 
thermodynamic efficiency of the engine cycle, thus increasing the range of the propulsion 
system. One of the critical hurdles to achieving detonation-based propulsion systems is 
creating an efficient deflagration-to-detonation process. Understanding the initiation 
mechanisms for an efficient detonation process is vital in overcoming the current 
technical challenges associated with implementing this type of propulsion. Current 
research at the Naval Postgraduate School (NPS) along with several other institutions are 
attempting to overcome the technical hurdles that are currently inhibiting the 
implementation of this advanced propulsion system. 

The primary advantage offered by pulse detonation engines (PDEs) is the high 
efficiency in a nearly constant-volume combustion process. Under Mach 5, specific 
impulses of PDEs have been theorized to be greater than most supersonic air-breathing 
engines currently in use, with the exception of advanced turbofan/turbojet engines. 
Eigure 1 shows various propulsion concepts and their specific impulse over a range of 
flight Mach numbers. 
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Figure 1. Comparison of High-speed Propulsion Teehnologies (From [1]). 

Figure 1 reveals that turbojet engines exhibit superior performance at Mach 
numbers below Mach 3, since above this value the complexity of the required advanced 
materials needed for these turbines minimizes the advantages of these propulsion 
systems. Above Mach 3, Ramjet and scramjet engines approach the performance level of 
turbojet engines. These engines do not have the ability to operate at low Mach numbers 
since the inlet air is compressed by the diffuser, which requires a booster to attain the 
initial flight speed necessary for engine start. 

Detonation-based propulsion systems, such as PDFs, have a greater operating 
range covering the subsonic region through high Mach numbers. The PDE combines 
high performance, efficiency, and simplicity thus making it a viable mode of propulsion 
for tactical missile applications. Other concepts incorporating detonation-based 
propulsion are a hybrid PDE gas turbine system and a combined cycle PDE for single 
stage to orbit launch vehicles. 
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In order to capitalize on the benefits of detonation-based propulsion, efficient 
initiation of the detonation in a short distance is critical. However, using obstacles to 
initiate detonation can result in total pressure losses, which results in a loss of efficiency. 
Comprehending the initiation mechanisms for detonation will allow designs that are more 
efficient and also maximize the frequency of repeatability. This research examines the 
initiation mechanisms through experimental observations and computer simulations for 
several obstacle configurations that are known to produce detonations in order to generate 
guidelines for the use of specific obstacle configuration. 
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II. BACKGROUND 


A. COMBUSTION PROCESSES 

Detonation is one of the three primary modes of eombustion. It is different from 
the other modes of eombustion based on the speed of propagation, the shoek wave 
ereated, and the effeets of this shook on the continuing process. A detonation event can 
be produced almost instantaneously from several options such as the use of an explosive 
charge, high energy ignition, or shock focusing if sufficient conditions exist. A 
conventional burning process can also be manipulated to transform into a detonation. 

1. Deflagration 

The most common form of combustion is deflagration. This method can be 
approximated by a constant-pressure combustion process; this form of combustion takes 
place in turbojet engines and also closely approximates diesel engines. In a perfect 
deflagration event, a combustion wave propagates through the mixture at subsonic 
speeds, with combustion and the resulting energy release occurring only at the flame 
front and leaving combustion products in its wake without a pressure rise. A more 
realistic model is that the deflagration wave speed is limited due to physical limitations 
that keep the velocities subsonic. The temperature, pressure, and structure of the 
advancing flame front have a robust affect on the chemical reaction rates of the 
substances involved. The chemical reaction rates govern the energy release rate that will 
further limit the flame speed, since combustion can only occur when the flame front 
makes contact with the reactants. Introducing turbulence into the event will increase the 
total flame area and subsequently increase the rate of consumption of the reactants. 
Deflagration results in a relatively large release of entropy, resulting in lower 
thermodynamic efficiencies [2]. 

2. Detonation 

Detonation occurs when a supersonic shock wave propagates through a reactive 
mixture and compresses it, which results in a prompt increase in temperature, pressure. 
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density, thus intensifying the reaction rate. An extremely energetic mixture develops 
directly behind the shock wave and promptly combusts in a violent, exothermic reaction 
that supports the previous shock wave and the two become dependent on each other. The 
compression and coupling properties in a detonation event result in a much larger 
enthalpy gain than from a deflagration event. 

3. Deflagration versus Detonation 

Comparing the attributes between deflagration and detonation emphasizes the 
differences between the two. A one-dimensional (ID) model of a combustion wave is 
shown in Figure 2 and reveals the basic properties of a combustion process. In the 
combustion wave reference frame, unburned reactants with velocity ui are located to the 
left of the stationary combustion wave and burned reactants with velocity U 2 are to the 
right. 


Stationary Combustion Wave 



^777777777777777777/777777Z777777777777/7777T/ 


Figure 2. Stationary One-Dimensional Combustion Wave Model (From [2]). 


Based on Figure 2, Table 1 lists the differences between deflagration and 
detonation for velocity (u), density {pj, temperature (7), and pressure (p) before and after 
the combustion wave. When compared with deflagration, the difference in pressure and 
temperature following the detonation wave are significant; this difference and the lower 
overall increase in entropy make detonation much more efficient than other types of 
combustion. 
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Table 1. Comparison of Detonation and Deflagration Characteristies (From [2]). 



Detonation 

Deflagration 

U\lc\ 

5-10 

0.0001-0.03 

Ullu\ 

0.4-0.7 (deceleration) 

4-16 

Pllp\ 

13-55 (compression) 

0.98-0.976 (slight expansion) 

TilT, 

8-21 (heat addition) 

4-16 (heat addition) 

pilpi 

1.4-2.6 

0.06-0.25 


B. DETONATION THERMODYNAMICS 

In order to understand the purpose of this researeh, a review of thermodynamics 
and detonations is necessary. The simple one-dimensional model of Figure 2 can be 
characterized thermodynamically through the Hugoniot curve. The assumptions applied 
to this elementary model include the perfect gas law as well as the conservation of mass, 
momentum, and energy. 


Perfect Gas Law: 

p = pRT 

(1) 

Conservation of Mass: 

P\^\ ~ Pi^i 

(2) 

Conservation of Momentum: 

PlUl = Pl+ Piul 

(3) 

Conservation of Energy: 

^ 1 2 ^ 1 2 
C/ 1 +-W 1 +^=c/2+-W2 

(4) 


where q is the specific heat added to the system and is the specific 
pressure, equivalent to: 


Where gamma is the ratio of specific heats: 


heat at constant 

(5) 


( 6 ) 
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For a constant area problem, Equation (2) illustrates that mass flow rate (m ) must 
be constant and when combined with Equation (3) yields the Raleigh-Eine relation. The 
slope of the Raleigh-Eine indicates the velocity of the detonation wave. 


Raleigh-Eine Relation; 


P]U\ 


Pi-Pi 

1 1 


= m 


P\ Pi 


(V) 


The Hugoniot Relation may be obtained by manipulating the Raleigh-Eine 
relation of Equation (6) through substitution of Equations (1), (3), (4), and (5) and 
solving for this yields a form that cancels the velocity conditions ul and u2, relating 
the heat release per unit mass. 


Hugoniot Relation; 


Y 


r-\ 


( 

El. 

V Pi 


\ 

El 

Ay 


1 


iPi-Pi) 


1 1 
- 1 - 

V A Pi J 


= q (8) 


Equation (8) describes the Hugoniot Curve of Eigure 3, which, for any 
combustion mixture with a given q and initial conditions p\, p\, displays all 
mathematically possible combinations of pi and Mpi. 



8 













Figure 3 states that there are five possible regions on the Hugoniot eurve, eaeh 
one a theoretieal eombustion eondition. These five regions are ereated by a eonstant 
pressure line, a eonstant speeifie volume line, and tangent lines drawn from the origin of 
the eurve defined by the initial eonditions p\ and p\ (labeled point A on Figure 3). The 
tangent lines interseet the eurve at the upper and lower Chapman-Jouguet (C-J) points, U 
and L. All of the regions, while mathematieally possible, are not observed; Region V is 
physieally impossible as it requires p 2 > pi and \/p 2 > Hpi, whieh would result in an 
imaginary veloeity for ui in Equation (7). 

Regions I and II are transient in nature. Region I is eonsidered the strong 
detonation region sinee p 2 > Pu, and requires an overdriven shoek. Region II is ealled the 
weak detonation region and requires exeeptionally fast ehemieal reaetions. Region IV is 
the strong deflagration region and is not a possible steady-state solution. Therefore, the 
only possible regions for steady state solutions are Region III, weak deflagration and the 
Upper C-J point, U, the solution point for detonations [2]. 


The veloeity of the burned reaetants ean be found by differentiating and re¬ 
arranging the Hugoniot Relation, resulting in the slope of the eurve. Equation (9): 


Slope of the Hugoniot Curve: 


dp2 


iPi-Pi)- 


^ 2/ ^ 

vr-iy 


Pi 


d{\/p,) f 2r ^ 
vr-ly 


J_ 

A 


V A Ay 


And equating it with the slope at the C-J points. Equation (9): 

dp 2 _ P2-P1 


Slope at C-J Points: 


d{l/P 2 ) 1/P 2 -I/Pi 


= K 


det 


Einally eombined with the Rayleigh-Line Relation of Equation (7): 
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Both detonation and deflagration results in the burned reaetants moving at a 
velocity equal to the local speed of sound and away from the wave. Directly behind the 
shockwave is a reaction zone that is characterized by subsonic flow, which allows the 
chemistry 

C. DETONATION STRUCTURE 

After the formation of a detonation wave, the event is maintained due to the strong 
coupling of the detonation shock wave and the combustion event that occurs just behind 
the shock wave. Zeldovich, von Neumann, and Doring developed this theory first for 
detonation and created a one-dimensional model that demonstrates this event [2]. Figure 
4 is a ZND model representing a one-dimensional detonation wave. 



Figure 4. ZND Detonation Wave Profile (From [2]). 


This ZND model consists of a leading supersonic shock wave, which causes a 
rapid increase in temperature and pressure in the reactants. After the shock wave is an 
induction zone; this reaction is traveling at subsonic velocities thus allowing the 
chemistry to support the initial shock by creating further sonic disturbances. A highly 
energetic and efficient chemical reaction takes place in the reaction zone, which sustains 
the shock wave. 
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The one-dimensional model is sufficient for explaining the critical characteristics 
of a detonation but the actual event is much more complex and is shown in Figure 5. 
Figure 5 uses a soot foil technique that is able to image the shock waves. The leading 
shock wave is actually comprised of many small, divergent shock waves that are 
followed by perpendicular shocks interacting with each other as well as the leading shock 
wave. These shock wave interactions form triple points, which are visible in Figure 5. 



Figure 5. Image of Three-Dimension Nature of a Detonation Wave (From [2]). 

D. THERMODYNAMIC ADVANTAGES OF DETONATIONS 

The useful energy lost in a thermodynamic process corresponds to an increase in 
entropy(s); the lower the rise in entropy due to combustion, the more energy available 
that can be extracted into useful work, thus a more thermodynamically efficient process. 
Relative to other combustion processes, detonation shows evidence of a very high 
thermodynamic efficiency. 

Figure 6 shows the relative values of entropy for the different regions of the 
Hugoniot curve of Figure 3. From Figure 6, it can be seen that the Lower C-J point 
corresponds to the maximum value of entropy which is described as strong deflagration. 
The minimum entropy value happens at the Upper C-J point, where detonation is the 
dominant process. Based on this observed difference in entropy for the various methods 
of combustion, detonation is the most efficient combustion process [2]. 
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Figure 6. Entropy Distribution on the Hugoniot Curve (From [2]). 

E. DEFLAGRATION-TO-DETONATION TRANSITION 

Detonation can be created by several initiation methods. Some of these methods 
include high-energy ignition, shock focusing, explosive charges or deflagration-to- 
detonation transition (DDT) [2], In order for detonation to become a viable option for 
propulsion, the detonation event must be repeatable and at a high enough frequency in 
which an effective amount of thrust is created. Initiating the detonation using DDT does 
not limit the frequency of repeatability and is also a relatively simple approach. 

Initially the deflagration process is started by a reactive mixture igniting and 
forming a small kernel of combustion. The resulting flame front then expands as it moves 
down the combustor and produces pressure waves ahead of the flame front. Ultimately, 
the multiple compression waves will combine into a single shock wave which results in 
the flame front breaking up due to the turbulence. Breaking up of the flame front causes 
an increase of the flame surface area thus increasing the reaction rates and energy release 
until an “explosion within an explosion” occurs in the reaction zone (Figure 7). Two 
resulting shock waves result: the super-detonation wave traveling forward into the 
unbumed gases and a retonation wave moving toward the rear of the combustor. In 
addition to these shock waves, a spherical shock wave develops as well as transverse 
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waves created from the oscillations between the two larger shock waves (Figure 8). A 
final steady detonation wave is created by the interactions of these shock waves as they 
travel down the combustor [2], 



Spherical 

Shock 



in Explosion" 

Figure 7. DDT “Explosion within an Explosion” (Erom [2]). 
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Eigure 8. DDT transverse Waves (Erom [2]). 
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Wall effects and the autoturbulization which results in high-intensity turbulence 
can cause the transition to detonation as long as the tube is of sufficient length. This 
autoturbulization causes additional wrinkling of the flame, which increases the flame 
surface area and leads to detonation. 

In order to shorten the distance required for DDT, turbulence can be intentionally 
created by placing obstacles in the flow field as shown in Figure 9. Obstacles in the flow 
field assist in the formation of shock-focused “hot spots” that lead to self-ignition of the 
fuel ahead of the flame front which results in the acceleration of the reaction zone. 
Together these phenomena considerably raise the reaction rates, which result in an 
increased compression wave strength. Obstacles currently being evaluated for their 
benefits include ramps, screens, tabs, and Shchelkin spirals. 
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Figure 9. Sequence of DDT Acceleration in Tube with Obstacles (From [3]). 
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F. 


PULSE DETONATION ENGINES. 


In order to produce thrust sufficient for propulsion, the detonation events must 
occur at a high enough frequency. This concept forms the basis of operation of the pulse 
detonation engine. Figure 10 demonstrates a single cycle of a simplified standard pulse 
detonation engine. This model assumes that a valve gas delivery system is used to direct 
reactants during the appropriate time. 



1. 

Reactants injected 

5. 

Detonation wave exits 

2. 

Tube fills 

6. 

Rarefaction waves 

3. 

Mixture ignited 

7. 

Products exhausted 

4. 

Detonation wave forms 

8. 

Cycle repeats 


Figure 10. Ideal PDF Operation Cycle (From [4]). 


Thrust is generated only during the detonation phases of this cycle; in order to 


produce a near-constant thrust, the overall cycle time must be held to a minimum. The 
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fill, blow-down, and purge cycle times are all extremely reliant on the length of the 
chamber. The obstacles that are used to enhance DDT can have a negative impact on the 
process by increasing the pressure loss through the combustor. Therefore, the ideal design 
for a PDE is one that minimizes the pressure losses from obstacles and the ignition delay 
due to filling requirements. 
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III. DESIGN/EXPERIMENTAL SETUP 


The experimental testing was conducted at the Naval Postgraduate School, Rocket 
Propulsion Laboratory, Test Cell #1. The test apparatus was comprised of a combustor 
section, shock formation tube, test section, optical viewing section, and exhaust tube 
(Figure 11). Gaseous ethylene (C 2 H 4 ) and air (O 2 + 3.76 N 2 ) were used as reactants and 
were pre-mixed prior to injection into the combustion section. The ignition of the 
ethylene/air mixture was initiated by a high capacitance igniter mounted in the head 
flange of the combustion section. The instrumentation used during testing included high¬ 
speed pressure data acquisition integrated with a high-speed digital camera that captured 
the images in the optical test section. The entire test assembly was supported by a 
Newport Research Corporation optical table. The actuation and instrumentation for the 
test cell was accomplished remotely in the control room via a National Instruments GUI. 
Operation of the camera was also controlled from the control room but from a separate 
computer. 



BANK 


Figure 11. Diagram of the Experimental Assembly. 
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Figure 12. Test Cell #1 Experimental Setup. 


A. DESCRIPTION OF THE EXPERIMENTAL ASSEMBLY 

I. Combustor 

The purpose of the eombustor seetion was to generate a fully developed 
deflagration wave and the assoeiated shoek. The eombustor seetion was made of a 
112.71 mm diameter stainless steel tube 340 mm in length (Figure 13). The inlet for the 
air/fuel mixture was located 25.4 mm from the combustor head and entered through two 
25.4 mm diameter tubes welded to the steel tube 180 degrees apart. Located on the head 
of the combustor was an igniter flange for the high capacitance igniter and it was shielded 
by a cone shaped shroud as shown in Figure 14. The purpose of the shroud was to 
protect the nascent flame from being extinguished during high flow rate testing. 
Attached to the downstream end of the combustion section is the shock formation tube. 
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Figure 13. Forward and Side View of the Combustion Section. 
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Figure 14. View of Igniter Flange with Shroud and Igniter Installed. 
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2 . 


Shock Formation Tube 


In order for the combustion shock to become fully developed and separated from 
the combustion wave, a shock formation tube of sufficient length was necessary. The 
shock formation tube enabled the formation of a combustion-driven normal shock wave 
prior to entering the test section. The shock formation tube was a 123.5 cm long, 13.0 cm 
square stainless steel tube. Michael A. Fludovich constructed this section for use in 
earlier thesis work [5]. To prevent deformation of this section during testing, four ribs 
made from 2.54 cm aluminum plate spaced along the shock formation tube. The test 
section was fastened to the aft end of the shock tube via fasteners and an o-ring seal. 

3. Test Section 

After the shock formation section was the test section which was also designed by 
Fludovich to allow a smooth transition from the larger 13.0 square shock formation tube 
to the 114 mm x 80 mm test section [5]. The test section was designed with several types 
of connections to allow various instrumentations. During this testing, two high-speed 
Kistler pressure transducers spaced approximately 10 cm axially along the centerline 
were installed as shown in Figure 15. High-speed Kistler pressure transducers were used 
to determine shock speed entering the test section. An optical sensor that was able to 
detect the flame front passage was also installed in this section, providing a trigger signal 
to the high-speed camera as shown in Figure 15. 
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Figure 15. Test Seetion Detail Showing Kistler Probes and the Optical Sensor. 



Figure 16. Test Section Drawing with Dimensions. 


4. Optical Test Section 

The optical test section was a previous nozzle inlet that was designed by 
Fludovich [5]. The nozzle inlet had been modified since the nozzle components had been 
removed and replaced with plates extending from the test section. The two windows in 
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the optical viewing section were 30.48 mm thick and were mounted flush with the 
sections walls, forming an imaging area 9 cm x 20.32 cm (Figure 16). The aft end of the 
optical test section had mounting brackets for attachment of the exhaust tube section. To 
facilitate additional testing, one of the windows was removed and replaced with an 
aluminum insert of the same dimensions. This aluminum insert allowed swept ramps to 
be mounted centerline on the vertical wall. This configuration allowed the existing 
imaging equipment to catch a top-down perspective of the detonation event as it occurred 
behind the obstacles. 

5. Ramp Mounting Plates 

Two ramp mounting plates were designed to allow mounting of the swept ramps 
both on the top and bottom of the test assembly. The mounting plates were 137.8 cm x 8 
cm and were tapered at the front end as shown in Figure 17. The two ramp mounting 
plates were mounted to the top and bottom of the test assembly via fasteners. The ramp 
mounting plates were mounted in four sections of the assembly: the shock formation tube 
section, test section, optical test section, and the exhaust tube. Both of the plates were 
constructed with an array of tapped holes which permitted multiple ramp configurations 
to be tested. The fastening holes were plugged when not in use to prevent any of the 
combustion products from filtering through the plate during testing. The same swept- 
ramp obstacles as used by Anderson [7] were used in this thesis. The details of the ramp 
geometry are given in Appendix A. With both ramp mounting plates installed, the 
resulting area for the optical test section was 76.2 mm x 114 mm (Figure 18). 
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Figure 17. Swept Ramp Removable Mounting Plate with 6 Rows of 2 Ramps. 
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Figure 18. End View of Combustor with Exhaust Tube Removed. 


6. Exhaust Tube 

A removable exhaust tube was installed downstream of the optical test section to 
isolate the test section from ambient flow disturbances and also direct the shock and 
combustion products away from the test table and instrumentation (Eigure 19). The 
exhaust section was a square tube which had two Kistler pressure transducers mounted 
centerline 10.16 cm apart along the right side of the exhaust tube. These sensors 
provided precise measurements of the shock speed after the test section which allowed a 
comparison to the upstream shock velocity measurement. The exhaust tube was inserted 
into the test cell exhaust tube, which minimized the pressure waves that could be directed 
back toward the instrumentation and cameras during testing. 
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Figure 19. Exhaust Tube and Test Cell Exhaust Tube. 


B. AIR AND FUEL DELIVERY 

High pressure (HP) air was provided to the test eell from the facility air system 
and ethylene was provided from a single bottle. The main control room computer 
controlled the supply pressures for both HP air and ethylene with the Tescom ER3000 
Version 2.0 software; this allowed the pressures of the reactants to be set remotely and 
independently of each other (Eigure 20). 
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Figure 20. Schematic Diagram of Air and Fuel Delivery System (From [8]). 


Air was supplied through PA-inch piping that was connected to a single inlet 
manifold and from there routed to the combustion section via flexible steel hoses. 
Ethylene gas was supplied via A-inch tubing that was connected to each of the air inlets 
prior to entering the test section; these are shown in Figure 21. 
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Figure 21. Fuel and Air Supply Delivery System. 


Air flow was initiated by the using the ER3000 regulator controller. Ethylene was 
injected into the combustor through a single Swagelock ball valve that was opened at the 
desired time for the desired time length. This ball valve was air actuated and controlled 
by the EabVlEW software via Crydom control solenoid switches located in an electronics 
cabinet located in the test cell (Eigure 22). 
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Figure 22. PXI-IOOOB Chassis (Upper Section) and Crydom Control Solenoid Switches 

(Lower Section). 

An orifice was located in both of the supply lines to provide an accurate method 
of controlling the mass flow rates for each gas. The choke diameters used during this 
testing were 0.236 inches for the air supply and 0.157 inches for the ethylene supply. A 
check valve located between the combustor and the ball valve was installed to prevent the 
possibility of backflow due to the combustion event (Figure 21). 

C. IGNITION SYSTEM 

The ignition system used in this test assembly was the Unison Vision-50 Variable 
Ignition System and is shown in Figure 23. This ignition system is a capacitive discharge 
type system that uses an aviation grade spark plug mounted in the head of the combustor 
section. Prior to testing the igniter was configured manually to provide a 2 Joule spark 
and was remotely triggered by the Lab VIEW software. 
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Figure 23. Unison Ignition System. 


D. INSTRUMENTATION 

The instrumentation used during testing controlled the initiation as well as data 
collection. Each of the four pressure signals were routed through one of four National 
Instruments (NI) 14-bit PXI-6II5 cards mounted in the NI PXI-IOOOB chassis, shown in 
the upper half of Figure 22. This chassis interfaced with the computers in the control 
room through the NI PXI-MXI-4 PXI Bridge and was capable of collecting either real 
time or high-speed buffered data. Optical data imaged from the high-speed camera was 
routed to a dedicated desktop computer in the control room using a fiber optic line. A 
schematic of the data acquisition system is shown in Figure 24. 
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1. Dynamic Pressure Transducers 

Before and after the test section were four high-frequency Kistler dynamic 
pressure transducers mounted in pairs in order to calculate the wave velocity based on the 
pressure transients. The pressure transducers were connected to Kistler Type 5010 Dual 
Mode amplifiers which then routed the data to two NI PXI-6115 data cards in the PXI- 
lOOOB for high-speed data collection. Once the testing sequence was initiated, the cards 
were configured to begin collecting data with the triggering of the capacitive discharge 
igniter and collected data at a rate of 500 kHz. The data was stored in the card buffer and 
was later saved to the computer located in the control room. The pressure data supplied 
by the Kistler transducers provided high resolution of shock passage. 
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2. Optical Sensor 

An optical sensor was inserted in the test section, 13.8 cm prior to the optical test 
section, and was used to activate the imaging equipment once it detected the passage of 
the flame front. With the flame detected, the optical sensor triggered the camera to begin 
colleetion of high-speed flame images. 

3. Ultra 17 High-Speed Imaging System 

An Intensified Ultra 17 CCD camera from DRS Data & Imaging Systems, Inc. 
was used to image the optical test section. The Ultra 17 was able to image and store up to 
17 high-speed frames per test run with a resolution of512x512 pixels per frame with a 
12 bit dynamie range. The Ultra 17 was equipped with an image intensifier which had 
the ability to provide a gain of up to 15,000 times in 100 discrete steps. The system was 
capable of imaging up to 157,000 frames per second [9]. 

The Ultra 17 eamera was controlled from the control room with a dedicated 
computer using a Hotlink fiber optic connection. Images from the test section were 
reflected off a mirror to the camera that was aligned parallel to the combustion tube as 
shown in Figure 25. The camera was triggered remotely from the signal generated from 
the optical sensor. A Nikon ED AF NIKKOR 80-200 mm 1:2.8 telephoto lens was 
attached to the camera allowing focused image of the optical test section (Figure 26). 
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Figure 25. Photo and Diagram Illustrating Camera Setup for Flame Imagery. 



Figure 26. Intensified Ultra 17 CCD High-speed Camera. 
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The Ultra computer program allowed the camera to be configured and armed 
remotely. The images were transferred to the control room computer running the Ultra 
software where the images were cataloged and stored. 

4. Princeton Instruments Camera 

Another camera utilized was a Princeton Instruments camera. The camera utilizes 
a charge-coupling device (CCD) array. This camera was triggered by an external optical 
sensor which then used a programmable gate pulse generator to input the appropriate 
delay. Mounted to the lens of the camera was a filter that allowed wavelengths of 425 nm 
to 435 nm to pass through. This wavelength correlates to chemiluminescence from the 
Carbon Hydrogen radicals (CH*) during a combustion event. 

5. Shimadzu Hyper-Vision 2 High-Speed Camera 

Approximately three quarters through testing the Ultra 17 malfunctioned and DRS 
Data & Imaging Systems loaned the lab a replacement. The replacement camera was a 
Shimadzu Hyper-Vision 2 high-speed camera (Figure 27). The camera was to image up 
to one million frames a second and store 99 frames. The camera speed was set to 125 
kfps allowing a 4 psec exposure time. A Nikon ED AF NIKKOR 80-200 mm 1:2.8 
telephoto lens was attached to the camera allowing focused image of the optical test 
section. The camera was triggered remotely from the signal generated from the optical 
sensor. 
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Figure 27. Shimadzu Hyper-vision 2 High-speed Camera. 


E. SOFTWARE DESCRIPTION AND FUNCTIONS 

The software used during this experimental testing included Cequel, Lab VIEW, 
and Ultra (Figure 28). These programs offered the ability to control all of the test 
instrumentation and equipment, but also provided the analysis tools. All applications 
were run on desktop PCs using Microsoft Windows. 


I. CEQUEL 


The chemical reaction characteristics were calculated using the CEQUEL 
Toolbox: Chemical Equilibrium in Excel Version 1.75 from Software and Engineering 
Associates, Inc. / Spreadsheet World, Inc. The CEQUEL plug-in was installed for use 
with Microsoft Excel 2003 and determined the theoretical combustion products and 
characteristics of different compounds for multiple reactions. 

The required inlet mass flow rates of air and ethylene supplied to the test section 
were determined from the desired equivalence ratio (cp). The equivalence ratio relates the 
fuel/oxidizer ratio of a given mixture to the ideal stoichiometric ratio of the mixture and 
is defined by Equation (12): 



^ J Stoichiometric 


( 12 ) 
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The stoichiometric ratio is a theoretical chemical reaction where complete 
combustion has taken place; therefore, no reactants remain after the event. In an actual 
combustion event, complete combustion will never occur since dissociation is always 
present. Another factor that prevents complete combustion is that complete mixing of the 
reactants does not occur, thus preventing the reactants in the mixture to react completely. 
The stoichiometric equation for an ethylene (C 2 H 4 ) and air (O 2 + 3 . 76 N 2 ) is shown in 
Equation (2): 


C 2//4 + 3(02 + 3.76A 2 ) ^2002 + 27/20 + 3 ( 3 . 76 ^ 2 ) 


(13) 


When a fuel/air ratio is fuel rich, unreacted fuel constituents will remain after 
combustion; this condition occurs when the value of the equivalence ratio (cp) is greater 
than one. When the equivalence ratio is less than one, the fuel/air ratio is fuel lean and 
unreacted oxygen remains after combustion. Regardless of the value of (p, nitrogen is 
present but is considered a non-reactive component. 


CEQUEL calculated the mass flow rates required to produce equivalence ratios 
ranging from 0.8 to 1.6. Inlet supply pressures for air and ethylene were determined by 
assuming a constant supply temperature for the choke diameters installed in each of the 
supply lines using the isentropic flow relationship of Equation (14) [10]; 


m = 


Choked Mass Elow Rate: 






(14) 


Kr = 


where; 
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(15) 


and Kr=l (choked flow). 

The air and ethylene pressures inputted into the ER3000 were calculated from 
Equation (14) with the mass flow rates previously determined from CEQUEE. 
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LabVIEW 


Actual testing was controlled remotely from the control room with LabVIEW 
Professional Development System Version 8.5 and a graphical Virtual Interface (VI) 
developed specifically for the experiment (Figure 28) installed on a desktop PC. When 
the LabVIEW software and the VI were placed in the run mode, the VI controlled the 
injection of fuel, igniter initiation, and collection/archiving of the test data. 



Figure 28. Screen Capture of LABVIEW VI Used During Testing. 

3. Ultra 17 Camera and Ultra Software 

The Ultra 17 high-speed camera was controlled from a dedieated desktop PC in 
the control room executing the Ultra Software, Version 1.1.19.1. The Ultra software is a 
Windows-based program issued by DRS Data & Imaging Systems, Inc. The software 
was used to program the imaging sequence as well as to arm the camera. Once armed the 
camera was triggered by a signal from the optical sensor located in the test section. 
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IV. EXPERIMENTAL RESULTS 


A. PURPOSE 

The purpose of this testing was to observe the eonditions before, during and after 
the detonation event. Various devices were used to record the actual event. The data 
recorded from the high-speed pressure transducers determined if a detonation occurred. 
The sequence of events surrounding a detonation event were recorded, analyzed and then 
compared to the results derived from the computer simulations in order to study the 
mechanisms responsible for the initiation of detonation. 

B. OBSTACLE FREE TESTING 

Initially testing was conducted with all obstacles removed from the test assembly 
in order to verify the set-up and operation of the imaging and data collection procedure. 
Multiple tests were conducted at various conditions to verily the test equipment was 
properly installed. The data recorded during all test included the pressure transducer 
values, fuel pressure, air pressure, the optical sensor output, and the images of the event. 
High-speed Kistler pressure transducers were used to record the pressure waves (shocks) 
produced by the combustion event due to the initial shock wave and then the actual 
combustion wave (Figure 29). 
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Figure 29. Fligh- Speed Pressure Traces for Obstacle Free Configuration. 

The combustor fill rate for this test had a characteristic Mach number value of 0.1 
with the fuel/air mixture at an equivalence ratio (cp) of 1.4. The pressure transducers 1 
and 2 recorded an initial pressure spike at 39.3 and 39.4 milliseconds, respectively, 
caused by the passing of the initial shock created by the combustion event (Figure 30). 
Knowing the distance between the probes and this time difference, a shock speed of 
558.3 m/s was calculated. The initial shock continued toward the exit of the assembly and 
registered on the pressure transducers 3 and 4 at 39.6 and 39.8 m/s respectively, which 
corresponded to a velocity of 561.3 m/s. The actual combustion wave follows the shock 
wave by approximately 0.75 milliseconds, which causes the second pressure spike in 
pressure transducers 1 and 2 at 39.95 and 40 milliseconds. The data recorded during the 
obstacle free testing reveals that the initial edge of the combustion wave follows the 
trailing edge of the initial shock wave by approximately 200 psec. 
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Figure 30. Pressure Traces of Figure 29 Magnified to Show Event Edges 

This approach to calculate the combustion and shock wave velocities were used in 
later tests with obstacles to initially determine if a detonation event occurred. 

C. SWEPT RAMP TESTING 

Two swept ramp configurations were tested. These arrangements were selected 
based on the success of previous testing conducted [7]. The first obstacle arrangement 
consisted of various rows of ramps located along the centerline of the ramp mounting 
plate and a mirrored ramp located on the opposing ramp mounting plate. The second 
configuration tested consisted of multiple rows consisting of 2 ramps mounted side by 
side. The testing for these arrangements consisted of 5 and 6 rows of ramps with cp 
ranging between 1.0 and 1.7. 
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Figure 31. High-Speed Pressure Traees for a 5 Rows of Centerline/Opposing Ramps. 

Figure 31 is the high-speed trace illustrating detonation in combustor for 5 rows 
of ramps mounted centerline and opposing. Pressure transducer #4 registered pressures in 
excess of 1300 psi. The combustor fill rate for this test had a characteristic Mach number 
value of 0.1 with the fuel/air mixture at an equivalence ratio (cp) of 1.23 and a shock 
speed of 2309 m/s. These measurements reveal that an over-driven detonation occurred. 
An over-driven detonation event is not a steady state condition but is a transient event 
since the detonation event will move downward along the Hugoniot curve (Figure 3) 
towards the Upper C-J point, where it becomes a stable detonation event. 

D. HIGH-SPEED IMAGERY 

I. Side View of Flame Imagery 

Once detonations were verified from the shock wave velocities, high-speed 
imaging was used to analyze the flame front interactions in order to better understand the 
initiation of detonations. Figure 32 shows selected frames of a detonation event. These 
images are from the same test that produced the high-speed traces of Figure 31. 
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Figure 32. High-speed Imagery of Flame Front Interaetions 5 Swept Ramps Located 

Centerline and Opposing. 
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From Figure 32, bright illuminations were visible behind the 4*’’ and 5* swept 
ramps as can be seen at 136 psec and 168 psec. These illuminations were located 
approximately 0.57 step heights behind the ramp (Figure 33). These illuminations 
appeared to be responsible for the detonation initiation observed shortly thereafter. 



Figure 33. Illumination Step Fleight. 


2. Top-Down View of Flame Imagery 

In order to determine if the detonations were occurring along the centerline of the 
ramps, imaging of the detonation from a top-down perspective were needed. To facilitate 
the additional imaging, one of the windows from the optical test section was removed and 
replaced with an aluminum insert of the same dimensions. This aluminum insert allowed 
swept ramps to be mounted centerline on the vertical sidewall opposite from the camera. 
This configuration allowed the existing imaging equipment to view a top-down 
perspective of the detonation event as it occurred behind the obstacles. The bottom edge 
of the optical window bisected the vertical ramps equally thus only the top half of the 
swept ramps was visible. 

To ensure detonation was achieved, 5 rows of 2 ramps were mounted upstream on 
one of the ramp mounting plates. Three additional ramps were then mounted to the 
vertical sidewall. With this configuration, a fuel/air mixture (q)=1.27) was ignited and 
detonation verified by wavespeed calculation. At time 108 psec and 132 psec, bright 
illuminations were present along the centerline and directly behind the 2"‘^ and 3'^‘^ 
vertically mounted ramps (Figure 34). 
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Figure 34. Top-down View of Flame Interacting with Obstacles 


3. CH* Imaging 

The bright illuminations visible behind the obstacles could have been caused by 
several scenarios, such as compression of the gases or a combustion event. In order to 
validate that these sites were due to the onset of detonation, additional evidence was 
needed. Another camera equipped with a filter was utilized in order to compare filtered 
versus unfiltered images. A Princeton Instruments camera was triggered from the same 
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optical sensor and imaged a single frame. Knowing the delay programmed into the 
Princeton Instrument’s camera, the same point in time was imaged from two different 
cameras and analyzed. Mounted to the lens of the Princeton Instrument’s camera was a 
filter that allowed wavelengths of 425-400 nm to pass through. This wavelength 
correlates to CH chemiluminescence from a combustion event and is an excellent marker 
for the onset of detonation. Testing was conducted with six ramps mounted centerline on 
both the top and bottom of the mounting plate. Numerous testing was required in order to 
capture the detonation event located behind the obstacles. Comparing the same event 
using the different cameras validated that the bright light visible was due to the molecular 
emission (CH) and was not a shock compression event (Figure 35). 



Figure 35. Capture of Detonation Behind 6th Ramp (Left); Corresponding Image with 

CH Filter Installed (Right) (q)=1.19) 


E. SUMMARY 

The purpose of this testing was not to create detonations using different 
configurations but to observe the flame front interactions in order to better understand the 
early development of detonations. Detonations were verified by the resulting shock 
speeds calculated from the high-speed pressure traces. Different images revealed that 
bright hotspots were formed along the centerline of the obstacles and approximately 0.5 
step heights behind the ramps. Using the appropriate filter, images were taken that 
indicated that the source of the bright illumination was from a combustion event and not 
from other sources such as shock heating. To better understand the initiation mechanisms, 
the detonation images need to be compared to computer simulations of similar 
conditions. 
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V. COMPUTER SIMULATION 


Computational Fluid Dynamic (CFD) simulations were used to explore the flow 
field/chemistry interactions responsible for DDT. The goal of this portion of the researeh 
was to better understand the mechanisms required for detonation while using various 
swept ramp configurations. 

Due to the complex chemistry, gas dynamics, and the varying pressure and 
temperatures oecurring during the event, the exact details of the detonation are not 
entirely understood. This complexity has had the effect of hindering the development of 
aceurate eomputer models and simulations that would assist in this investigation. The 
CFD simulations presented in this paper are to mainly illustrate the effect the of the flow 
field around the obstaeles with and without combustion present. 

The ramp configurations selected for examination required massive amounts of 
proeessor time and memory. A single 3D model typieally had more than 9 million 
elements; due to the large number of input variables and the desired output parameters, 
the computing power of a single maehine was insufficient. Due to this a high- 
performance computing cluster eonsisting of 32 individual quad-eore maehines running 
the ANSYS software package CFX was used to deerease the solution time. The required 
time for solution varied greatly depending on the type of solution, steady state or 
transient type, and whether chemistry was being simulated in the model. 

A. MODELING SOFTWARE 

Two different software packages were used in the eomputer simulation proeess. 
SolidWorks 2008 was utilized and is developed by Dassault Systemes SolidWorks 
Corporation. The second software package employed was ANSYS CFX and was 
designed by ANSYS, Ine. Both of these software paekages run on a Windows operating 
system; additionally, CFX ran on the cluster that was a Linux operating system. 
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1 . 


SolidWorks 


SolidWorks 2008 was used to draw 3D combustor models used in the computer 
simulation. SolidWorks is a 3D mechanical CAD software package that is a parasolid 
based modeler and incorporates a parametric based approach used to create models. 

The combustor and swept ramp obstacles were drawn to accurately simulate the 
actual conditions of the experimental test assembly. Various combustor lengths were used 
depending on the ramp configuration in order to produce valuable simulations on a 
practical time scale. Based on successful experimental results, two different ramp 
configurations were selected for further analysis using computer simulation. 

The challenge in modeling the combustor was visualizing the volume available 
for the reactants to occupy. In order to create this type of geometric model a drawing of 
the internal volume of the combustor was drawn (Figure 36) as well as the desired swept 
ramp obstacle configuration (Figure 37). 



Figure 36. SolidWorks Drawing of the Internal Volume of the Combustor. 
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Figure 37. SolidWorks Drawing of 6 Swept Ramp Stages Located Centerline and on the 

Top and Bottom of the Combustor. 

With these two part drawings, a geometric volume comparison was conducted in 
SolidWorks which resulted in the total ramp volume removed from the internal volume 
of the combustor (Figure 38). 



Figure 38. Volume Comparison Completed from SolidWorks Used in Simulation. 
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2 . 


ANSYS CFX 


ANSYS CFX is a high performance, general purpose Computational Fluid 
Dynamics (CFD) program that is used to solve various fluid flow problems. ANSYS 
CFX uses a solver that can be highly parallelized and incorporates a flexible user 
environment that employs intuitive GUIs. 

After the combustor volume was completed in SolidWorks, it was exported into 
ANSYS CFX. Once imported into CFX-Pre, a volume mesh was generated based on user 
inputs which control the quality and size of the resulting mesh (Figure 39). 



Figure 39. Example of Volume Mesh Generated in CFX-Pre. 

Once an adequate volume mesh was created, CFX-Pre was then used to define the 
boundary and initial conditions of the flows to be evaluated. All the simulations defined 
the forward end as the inlet and the aft end as an outlet; all other surfaces were treated as 
a wall. With this defined simulation, CFX computed a solution of the entire test volume 
at a time interval defined by the user. CFX-Post was used to view the results in sequence 
and to conduct post-processing and analysis on the results. Analysis was conducted by 
creating various planes through the combustor and defining parameters of interest to be 
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evaluated. These parameters included but were not limited to velocity, pressure, 
temperature, and turbulence kinetic energy. Figure 40 shows an example of an XY plane 
cut through the Z-axis of a solution. 



Figure 40. Example of a Post-Process Temperature Analysis of an XY Plane Slicing 

Through the Origin. 

B. SIMULATION RESULTS 

I. Swept Ramp Configurations 

CFX simulations of several ramp configurations were conducted to obtain steady- 
state results. The configurations selected were ones that detonations were observed 
during experimental testing. The configurations considered included five or six swept 
ramp rows constructed into two different arrangements. The first arrangement consisted 
of either five or six rows of two (Figure 41) while the second layout was either five or six 
rows of a single ramp located centerline with an opposing ramp located on the opposite 
face (Figure 42). The obstacles used had a maximum height of 19.43 mm and a sweep 
angle of approximately 63 degrees. The exact dimensions of the obstacles used are found 
in Appendix A. The typical initial conditions in the test volume and at the outlet were set 
for approximately sea level conditions to include a static pressure of 1 atm and 
temperature of 450 K. Air modeled as an ideal gas was simulated to flow into the inlet 
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with an initial mass flow rate of 0.3125 kg/s. The wall was simulated as an adiabatie 
surface with no slip conditions established for flow. A summary of the CFX steady-state 
simulations and the parameters used for the initial and boundary conditions are included 
in Table 2 and Table 3. 



Figure 41. Cutaway of Simulation Configuration of 6 Rows of Ramps 

Centerline/Opposing. 
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Table 2. Steady-State CFX Configurations with Resulting Volume Mesh Information. 


Configuration 

Ramps 
per Row 

No of 
Rows 

No of 
Nodes 

No of 
Elements 

5 Rows of 2 

2 

5 

2475988 

8756895 

6 Rows of 2 

2 

6 

2454757 

8709301 

5 Rows 

Centerline/Opposing 

1 

5 

2605670 

8961295 

6 Rows 

Centerline/Opposing 

1 

6 

2621557 

8974794 


Table 3. Steady-State CFX Initial and Boundary Conditions. 


Location 

Type 

Conditions 

Inlet 

Inlet 

Elow Direetion: Normal to Boundary Condition 

Elow Regime: Subsonic 

Heat Transfer: Static Temperature 

Static Temperature: 450[K] 

Mass Elow Rate: 0.3125 [kg/ s] 

Mass And Momentum: Mass Elow Rate 

Eddy Eength Scale: 0.02 [m] 

Eractional Intensity: 0.01 

Turbulence: Intensity and Eength Scale 

Outlet 

Outlet 

Elow Regime: Subsonic 

Mass And Momentum: Average Static Pressure 

Pressure Profile Blend: 0.05 

Relative Pressure: 0 [Pa] 

Pressure Averaging: Average Over Whole Outlet 

Wall 

Wall 

Heat Transfer: Adiabatic 

Mass And Momentum: No Slip Wall 

Wall Roughness: Smooth Wall 


2. Steady-State Results 

ANSYS CFX was used to determine the steady-state results for the eonfigurations 
listed in Table 2. Post-proeessing on CFX allowed multiple planes to be inserted at 
various loeations of the eombustor displaying a wide variety of parameters. Turbulenee 
Kinetie Energy (TKE) as well as shear strain rate were used for insight in understanding 
the eauses of detonations sinee areas of high TKE and shear strain rate were eorrelated to 
very reaetive loeations observed during experimental testing. TKE is the average kinetie 
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energy of the fluetuating veloeity eomponents in turbulent flow. The shear strain rate of a 
fluid is defined as the rate of deerease of the angle formed by two perpendieular lines in 
the fluid. Areas of high shear strain rate aet as a viscous boundary and thus have the 
ability to direct/divert flow. Figures 43-47 were created from the steady state simulation 
with 6 stations of inline and opposing ramps (Figure 41). 

Turbulence Kinetic Energy (KExy) 

>*2 s '- 2 ] 




Figure 43. Centerline Turbulence Kinetic Energy of 6 Ramp Rows Inline/Opposing. 



Figure 44. Centerline Turbulence Kinetic Energy in the 4th-6th Stages. 
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Figure 45. Centerline Temperature Plot with Velocity Vectors of 6 Ramp Rows 

Inline/Opposing. 

Figure 43 and 44 are centerline plots of the TKE in the combustor. Theses plots 
reveal an increase of TKE directly behind the swept ramps which grows in intensity 
down the combustor. Figure 45 is a centerline plot of temperature, which shows an 
increase in temperature behind each ramp, also growing in intensity down the combustor. 





Figure 46. Shear Strain Energy with Velocity Vectors (Located at Base of 5th Ramp). 
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Figure 47. Shear Strain Energy with Velocity Vectors (Located 2.54 cm Downstream of 

5 th Ramp). 

Figure 46 and 47 are plots of shear strain rate at different locations in the 
combustor. Both figures reveal a high shear strain rates directly behind the ramps, caused 
by the symmetrical vortexes that are created by the flow field rolling off each edge of the 
swept ramp. 

Figures 48 through 50 are plots of the steady state simulation of 6 rows of 2 ramps 
side by side (Figure 42). Figure 48 shows the centerline TKE down the length of the 
combustor; this plot reveals an increase in TKE behind each ramp station, which also 
increases in intensity down the combustor. 
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Figure 48. Centerline Turbulence Kinetic Energy of 6 Rows of 2 Ramps. 

Figures 49 and 50 plot the TKE at various locations; both reveal areas of high 
TKE behind the ramps and each slanted toward the outboard of the combustors. These 
plots show the resulting flow flelds due to the steady state flow interactions between the 
two ramps and the center channel. The resulting flow flelds due to the ramps were 
symmetrical along each station. All of the steady state simulations indicated areas of high 
TKE and temperature behind each obstacle, increasing in intensity down the combustor. 



Figure 49. Turbulence Kinetic Energy with Velocity Vectors @ Halfway Up 2nd Set of 

Ramps. 
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Figure 50. Turbulence Kinetic Energy with Velocity Vectors Midway between 5th and 

6th Ramps. 

3. Transient Results with Chemistry Simulated 

CFX simulations of several ramp configurations were conducted to obtain 
transient chemistry results. The transient simulation was initialized with the Ethylene/Air 
mixture model including all the products of combustion. Nitrogen gas was suppressed, 
which restricted the combustion to a modified ethylene/oxygen mixture. Combustion was 
modeled using finite rate chemistry including the eddy dissipation model and turbulence 
was computed with the standard k-epsilon model. The reference pressure and temperature 
were set to standard conditions. The reduced ethylene/air chemistry model used was 
based on that used by Foreman Williams of the University of California San Diego 
(UCSD) [11]. Initial inlet boundary total pressure was 2.7 atmospheres and 2,800K static 
temperature with the products of combustion being introduced at 600 m/s. The outlet 
boundary was set to subsonic conditions at 3.0 atmospheres and the composition was a 
mixture of ethylene and oxygen. Initially the combustor was filled with an ethylene/ 
oxygen mixture, once the simulation was started only the designated products were 
injected into the inlet. By injecting only the products, an ideal combustion event was 
created with only the designated products remaining in the wake of the combustion wave. 
An Arrhenius reaction rate with adequate activation energy was assigned to ensure a 


58 


















proper rate of reaetion took place. The transient simulation used an adaptive time step 
using RMS Courant number. To ensure adequate convergence, a value of four for the 
Courant number was chosen for these simulations. A maximum of 75 iterations per time 
step was selected to allow the simulation to converge. 

Based on the complexity of the simulation, the entire obstacle model could not be 
simulated. Since the entire obstacle configuration could not be analyzed, detonations 
were not generated in the simulation, but the fuel/air mixture flow fields could be 
examined at locations where the experimental results revealed strong reaction sites. In 
order to get a representation of the mechanisms required for initiating detonation and to 
simplify the simulation, two configurations were analyzed. A summary of the CFX 
transient simulations and the parameters used for the initial and boundary conditions are 
included in Table 4 and Appendix B. The first transient simulation consisted of one 
station of 2 ramps mounted side by side mounted on the bottom of the test area (Figure 
51). 



Figure 51. Transient Chemistry Simulation Model with 1 Station of Side by Side Ramps. 
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Table 4. Transient Chemistry CFX Configurations with Resulting Volume Mesh 

Information. 


Configuration 

Ramps 
per Row 

No of 
Rows 

No of 
Nodes 

No of 
Elements 

Single Station; 1 Row 
of 2 

2 

1 

2475988 

8756895 

Two Station; 2 Ramps 
Inline/Opposing 

1 

2 

2454757 

8709301 


Figure 52 represents the position of the flame front at various times as it passed 
over the middle of the ramp. The first time step behind the ramp revealed the absence of 
fuel since the flame front had been accelerated behind the ramps. Figure 53 are plots of 
the shear strain rate downstream of the ramp station and indicate the flow path. Figure 53 
revealed that a strong vortex was created due to the flow field rolling off the outboard 
edge of the ramps. A strong vortex was also created down the centerline of the combustor 
ramp located between the ramps; the interaction between these vortexes occurred directly 
behind the ramps as indicated by an area of high TKE which corresponded to observed 
regions of very rapid reactions (explosions). Figure 54 is a 3D representation of the 
vortexes created from the combustion event and from the interactions with the side by 
side ramp obstacles. 
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Figure 52. 



Ethylene Mass Fraetion Along the Centerline of a Ramp at Various Times. 
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Figure 53. 
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Figure 54. Ethylene Mass Fraetion Vortexes at Time 3.22E-4 see. 


After the results from this model were examined, a second arrangement was 
created in order to show the interactions between multiple stations of ramps. This second 
configuration consisted of 2 stages of inline ramps mounted along the centerline on the 
top and bottom of the test area (Eigure 55). Eor this second configuration, the number of 
elements in this model had to be reduced in order for the high-performance computing 
cluster to produce results. 
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Figure 55. Transient Chemistry Simulation Model with 2 Stations of Inline/Opposing 

Ramps. 

Figures 56 and 57 revealed a snapshot in time of the differing axial positions of 
the flame front. Figure 57 indicated that the flame front along the wall was traveling at a 
higher velocity that the flame front along the centerline. The flame front was accelerated 
behind the swept ramps due to the flow fields rolling off the edges of the swept ramps 
represented by the presence of combustion products present behind the ramps (Figure 
56). This difference in velocities explained the appearance in the images that the flame 
front had already passed the detonation site. A comparison was conducted between the 
simulation results and the images recorded from actual testing (Figure 58). This 
comparison revealed that flame front along the wall and the flame front behind the ramps 
were moving close to the same velocity, thus creating a 2D image that the flame front had 
passed the detonation site and then collapsed back toward it. 
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Figure 56. 


Figure 57. 
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Ethylene Mass Fraction along Centerline at Time 6.31E-4 sec. 



Ethylene Mass Eraction along the Wall at Time 6.31E-4 sec. 
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Figure 58. Comparison of Shape of Detonation Site with Aetual Image versus Ethylene 

Mass Fraetion along the Centerline. 

Figures 59 through 62 indicated that the rolling vortexes off each swept ramp was 
not symmetrical, but actually a dominant vortex off one side of each of the ramps caused 
a global rotation of the entire flow fields. The interactions between the flow field off the 
ramp edge and from the upstream ramps created areas of high TKE directly behind the 
ramps and correlated to observed detonation initiation locations in actual testing. 
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Figure 59. Ethylene Mass Fraction Vortexes at Time 6.3 lE-4 sec. 



Figure 60. Turbulence Kinetic Energy 0.4 cm behind 2nd stage at Time 6.3 lE-4 sec. 
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Figure 61. Shear Strain rate 0.4 em behind 2nd stage at Time 6.3 lE-4 sec. 
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VI. CONCLUSIONS AND FUTURE WORK 


A. CONCLUSIONS 

This thesis explored the fluidie/combustion meehanisms for swept ramp obstacles 
that are responsible for initiating a detonation event through experimental testing and 
computer simulation. The experimental testing revealed where the detonations occurred 
relative to the installed obstacles. The images revealed that the initiation sites for 
detonation occurred directly behind the ramps and within approximately one half of its 
step height. The orthogonal integrated view of the detonation event indicated that the 
flame front had passed the detonation site prior to the detonation event. This observation 
indicated the possibility of the flame front actually collapsing back into the detonation 
site and then the creation of the detonation. 

Computer simulations were created in an attempt to provide additional qualitative 
information about the observations from testing but more importantly create a simulation 
where multiple parameters could be analyzed at various locations in the flow fields that 
were not achievable through experimental testing. The initial simulations were without 
combustion and used the entire obstacle configuration. These simulations revealed the 
flow fields that were created from the interactions with the various obstacle 
configurations. The final type of simulations analyzed included the introduction of 
reduced chemistry for the ethylene/oxygen mixture. The simulations including chemistry 
indicated that the combustion process altered the flow fields at locations behind the swept 
ramps that were represented by areas of high TKE. After further analysis, the simulation 
revealed that the flame front did not collapse back toward the location of detonation. The 
results showed that the rolling flow field off of the swept ramps interacted with the global 
flow fields created sites that are favorable for detonation. 

B. FUTURE WORK 

The experimental results of this thesis show that areas of high TKE are precursors 
for detonation. Eurther refinement of the chemistry model and adjusting the maximum 
time steps will produce more precise results. A better method of developing an adequate 
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volume mesh is needed in order to faeilitate simulating the entire obstacle and also allow 
a more detailed investigation of the combustion events occurring after the ramps. 
Currently the mesh needed to produce results with adequate quality is greater than the 
capacity of the high-performance computing cluster. Future simulations could be 
conducted on different software thus allowing a more exact solution. 
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APPENDIX A: OBSTACLE DRAWINGS 


The drawing contained in this Appendix details the design of the swept ramp 
obstacle utilized in this research. 



1. Design Details for the Swept Ramp Obstacles [3]. 
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APPENDIX B: TRANSIENT SIMULATION INITIAL 

CONDITIONS 


Table 1. Domain Physics 


Domain - Default Domain 

Type 

Fluid 

Location 

B34 

Materials 

Ethylene Air Mixture 

Fluid Definition 

Material Library 

Morphology 

Continuous Fluid 

Settings 

Buoyancy Model 

Non Buoyant 

Domain Motion 

Stationary 

Reference Pressure 

l.OOOOe+00 [atm] 

Combustion Model 

Finite Rate Chemistry and Eddy 
Dissipation 

Eddy Dissipation Model 

Coefficient B 

5.0000e-01 

Component 

C2H4 

Option 

Automatic 

Component 

CO 

Option 

Automatic 

Component 

C02 

Option 

Automatic 

Component 

H 

Option 

Automatic 

Component 

H20 

Option 

Automatic 

Component 

N2 

Option 

Constraint 

Component 

02 

Option 

Automatic 

Component 

OH 
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Option 

Automatic 

Heat Transfer Model 

Total Energy 

Thermal Radiation Model 

P 1 

Spectral Model 

Gray 

Turbulence Model 

k epsilon 

Turbulent Wall Functions 

Scalable 


Table 2. Boundary Physics 


Domain 

Boundaries 

Default 

Boundary - 1 niet 

Domain 

Type 

INLET 


Location 

1 nIet 


Settings 


Component 

C2H4 


Mass Fraction 

O.OOOOe+00 


Option 

Mass Fraction 


Component 

CO 


Mass Fraction 

4.0873e-02 


Option 

Mass Fraction 


Component 

C02 


Mass Fraction 

1.3498e-01 


Option 

Mass Fraction 


Component 

H 


Mass Fraction 

O.OOOOe+00 


Option 

Mass Fraction 


Component 

H20 


Mass Fraction 

7.1727e-02 


Option 

Mass Fraction 


Component 

02 


Mass Fraction 

1.6936e-02 


Option 

Mass Fraction 


Component 

OH 


Mass Fraction 

7.7709e-03 


Option 

Mass Fraction 
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Flow Regime 

Mixed 


Blend Mach Number 
Type 

Normal Speed 


Fleat Transfer 

Static Temperature 


Static Temperature 

2.8000e+03 [K] 


Mass And Momentum 

Cartesian Velocity Components and 
Total Pressure 


Relative Total 

Pressure 

2.7000e+00 [atm] 


U 

6.0000e+02 [m s""-l] 


V 

O.OOOOe+00 [m s""-l] 


W 

O.OOOOe+00 [m s""-l] 


Thermal Radiation 

Local Temperature 


Turbulence 

Low Intensity and Eddy Viscosity 
Ratio 


Boundary - Outlet 


Type 

OPENING 


Location 

Outlet 


Settings 


Component 

C2H4 


Mass Fraction 

6.3744e-02 


Option 

Mass Fraction 


Component 

CO 


Mass Fraction 

O.OOOOe+00 


Option 

Mass Fraction 


Component 

C02 


Mass Fraction 

O.OOOOe+00 


Option 

Mass Fraction 


Component 

H 


Mass Fraction 

O.OOOOe+00 


Option 

Mass Fraction 


Component 

H20 


Mass Fraction 

O.OOOOe+00 


Option 

Mass Fraction 


Component 

02 
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Mass Fraction 

2.1815e-01 


Option 

Mass Fraction 


Component 

OH 


Mass Fraction 

O.OOOOe+00 


Option 

Mass Fraction 


Flow Regime 

Subsonic 


FI eat Transfer 

Opening Temperature 


Opening Temperature 

2.9800e+02 [K] 


Mass And Momentum 

Entrainment 


Relative Pressure 

3.0000e+00 [atm] 


Thermal Radiation 

Local Temperature 


Turbulence 

Zero Gradient 


Boundary - Wall 


Type 

WALL 


Location 

Wall 


Settings 


FI eat Transfer 

Adiabatic 


Mass And Momentum 

No Slip Wall 


Thermal Radiation 

Opaque 


Diffuse Fraction 

l.OOOOe+00 


Emissivity 

5.0000e-01 


Wall Roughness 

Rough Wall 


Sand Grain 

Roughness Fleight 

2.0000e+01 [micron] 
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APPENDIX C: CEQUEL CALCULATIONS AND TABLES 


CEQUEL was used to determine the mass flow rates necessary to achieve the 
various equivalence ratios during testing. Based on the required mass flow rates and the 
choke size installed in the piping, a pressure was determined that could be entered into 
the ER3000 for both air and ethylene. Both the air and fuel pressure gages were 
calibrated prior to testing to ensure that correct mass flow rate would be delivered in 
order to prevent any misinterpretation of data. 
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1. Table of Mass flow rates provided by CEQUEL. 
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Mtich ? 
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D5 
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AIR 

FUEL 

AIR 

FUEL 

AP 

Fuel coi 

CD 

0.8 

337.5 

32.7 

689.7 

80.0 

28.5 

108.5 

rn 

(N 

0.9 

335.5 

38.3 

685.6 

91.2 

31.9 

123.1 

o 

1.0 

333.4 

43.8 

681.6 

102.3 

35.0 

137.3 

"o 

1.1 

331.4 

49.3 

677.6 

113.3 

38.0 

151.3 

jz LD 

O 

1.2 

329.4 

54.7 

673.6 

124.1 

40.8 

164.9 

I_ o 
<!■ Q> 

1.3 

327.5 

60.1 

669.7 

134.8 

43.4 

178.2 

CO- o 

1.4 

325.6 

65.4 

665.8 

145.4 

45.9 

191.3 

o o 

1.5 

323.6 

70.6 

662.0 

155.9 

48.2 

204.1 


1.6 

321.8 

75.8 

658.2 

166.3 

50.4 

216.6 

5 ^ 
o 

1.7 

319.9 

80.9 

654.5 

176.5 

52.4 

228.9 


1.8 

318.0 

86.0 

650.8 

186.6 

54.3 

240.9 

o 

1.9 

316.2 

91.0 

647.1 

196.7 

56.0 

252.7 

< 

2.1 

312.6 

100.8 

639.9 

216.4 

59.1 

275.5 


2. Table of Air and Corrected Fuel Pressures Input to the ER3000. 


Fuel Pressure Transducer 



3. Fuel Pressure Transducer Calibration Curve for Fuel Pressure Correction. 
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Air Pressure Transducer 



4. Air Pressure Transducer Calibration Curve for Air Pressure Correction. 
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APPENDIX D: STANDARD OPERATING PROCEDURE 


Test Cell #1 Standard Operating Proeedures (S.O.P) 

Combustor Start Up 
(Modification Date 28 Oct 2009) 

INITIAL PREPARATIONS 


1. Notify all lab personnel of intention to make test eell 1 live. 

2. Turn ON eontrol console 

3. Turn ON warning lights 

4. Cell #1 EMERGENCY SHUTDOWN BUTTON (Control Room)- VERIFY IN 

5. Notify the Golf Course (x2167) (Only required if Hot Eire Test is conducted) 

6. Open Test Cell Door 

7. Igniter Control (Test Cell)- VERIFY OFF (Red Button Out) 

8. PXI- 1 OOOB Rack (Test Cell)- VERIFY ON 

9. Shop Air- VERIFY > 100 PSI 

10. Shop Air Valve (Test Cell)- VERIFY OPEN 

11. 115 VAC Control/Ceh #1 Switch (Control Room)-ON 

12. 28VDC Power Supply/Ceh #1 Switch (Control Room)-ON 

13. Open EABVIEW 


TESTING SET-UP 


1. Viewing Section-Disconnect from Pill Tube (Includes unbolting ski ramps) 

**Note: Viewing Section must be removed and slid aft in order to continue** 

2. Install the ramp obstacles 

3. Viewing Section-Attach to Pill Tube and bolt down ski ramps 

4. Kistler Amplifier Power-ENSURE OFF 

5. Kistler Eeads-CONNECT 

6. Exhaust Tube- VERIFY PROPER POSTION 

7. Notify ah personnel that gasses and TESCOM will be enabled. 

8. Test Cell #2 and #3 Node 1 Air Isolation Valve (Test Cell #2)- VERIFY 
CLOSED 

**NOTE: This valve maybe left open only if Test Cell #2 or #3 is configured for 

active testing** 


9. Turn on Optical Sensor. 

10. TRANSDUCER and TESCOM Power Switch (Test Cell #2)-ON 
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11. Set 0 (Zero) pressures on ER3000 (Control Room) for the following; 

a. Node 3 (Ethylene) 

b. Node 1 (Main Air) 

12. Main HP Air Jamesbury Valve (Outside Test Cell)-OPEN SLOWLY 

13. Power Strip (above PXI-IOOOB)-VERIFY ON 

14. Igniter Control Light (Red LED upper left on CRYDOM in PX-IOOOB Rack)- 

VERIFY OUT 

**DANGER: IF RED LIGHT IS ENERGIZED, MUST CLEAR USING LABVIEW 
BEFORE CONTINUING TO PREVENT PREMATURE IGNITION** 

15. Igniter Control (Test Cell)-PUSH RED BUTTON IN 

16. Igniter Control Startup Diagnostics-OBSERVE COMPLETION OF 
DIAGNOSTICS (Verify energy level setting reads 2.03 J) 

17. Main HP Air Isolation Valves (2) (Loeated in Test Cell)-OPEN SLOWLY 

**DANGER: OPEN VALVES SLOWLY TO PREVENT RAPID 
PRESSURIZTION OF DOWNSTREAM LINES** 

18. Node 3 (Ethylene) Shop Air Valve (Above Bottle in Bottle Room)-OPEN 

19. Ethylene Bottle Isolation Valve (On Bottle)-OPEN SLOWLY 

** VERIFY ADEQUATE ETHYLENE PRESSURE FOR TESTING ON 
DOWNSTREAM GAGE LOCATED IN BOTTLE ROOM** 

PRE-TESTING 

1. ULTRA Camera-TURN ON AND REMOVE COVER 

2. Start ULTRA Software 

3. Verify Image on ULTRA Software 

4. Verify Desired Trigger Types and Valve Durations 

5. Determine Desired Luel and Air Pressures (Mass Plow Choke Calibrations.xls) 

6. Set Required Pressures on ER3000 

a. Node 3 (Ethylene) 

b. Node I (Main Air) 

7. Insert Tape into VCR 

8. Switch Monitor to B (Test Cell #1) 

TESTING 

1. Clear All Test Cells and Verify with Head Count 

2. Plashing Red Light and Siren-ENERGIZE 

3. Verify Golf Course is CLEAR 

4. VCR-RECORD 

5. Camera-ARM 
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6. In LAB VIEW Enable Eacility Button-ON 

7. Test Cell #1 Emergeney Shutdown Button-TURN CLOCKWISE 

8. In LAB VIEW Start Button-CLICK TO START 

WHEN TESTING COMPLETE 

9. Set Node 1 Pressure to 0 (Zero) 

10. In LAB VIEW Turn Off Button-CLICK TO SECURE 

11. In LAB VIEW Enable Eaeility Button- VERIFY OFF 

12. Test Cell #1 Emergeney Shutdown Button-PUSH IN 

13. Siren-OFF 

14. VCR-Stop/Pause 

15. Save UETRA Image 

**NOTE: If Further Testing is Required, re-perform steps I-I5 of the Testing 

Section** 

POST-TESTING 


1. Set ER300 Node 3 (Ethylene) to 0 (Zero) 

2. Verify pressures on ER3000 are set to 0 (Zero) on the following: 

a. Node 3 (Ethylene) 
h. Node I (Main Air) 

3. Main HP Air Isolation Valves (2) (Eocated in Test Cell)-CLOSE 

4. Main HP Air Jamesbury Valve (Outside Test Cell)-CLOSE 

5. Ethylene Bottle Isolation Valve (On Bottle)-CLOSE 

6. Igniter Control (Test Cell)-PUSH RED BUTTON OUT 

**NOTE: If Further Testing will be accomplished with a different Ramp obstacle 
configuration, return to the TESTING SET-UP Section. If not continue to step 7.** 

7. ULTRA Camera-TURN OFF AND INSTALL COVER 

8. TRANSDUCER and TESCOM Power Switeh (Test Cell #2)-OFF 

9. Close Test Cell Door 

10. Node 3 (Ethylene) Shop Air Valve (Above Bottle in Bottle Room)-CLOSE 

11. Secure Bottle Room 

12. Exit out of ULTRA 

13. EXIT out of EABVIEW 

14. 28VDC Power Supply/Cell #1 Switch (Control Room)-OFF 
15.115 VAC Control/Cell #1 Switch (Control Room)-OFF 
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APPENDIX E: TEST SUMMARY 


Date 

C|J 

Wavespeed 

(m/s) 


Obstacle Configuration 

9/30/2009 



Bad Data 

6 rows of 2 

9/30/2009 

1.13 

2208.7 


6 rows of 2 

9/30/2009 



Bad Data 

6 rows of 2 

9/30/2009 



Bad Data 

6 rows of 2 

9/30/2009 

1.32 

2309.1 


6 rows of 2 

9/30/2009 

1.33 

4618.2 


6 rows of 2 

9/30/2009 



Bad Data 

6 rows of 2 

9/30/2009 

1.25 

1881.5 


6 rows of 2 

9/30/2009 

1.23 

423.3 


6 rows of 2 

9/30/2009 

1.33 

1539.4 


6 rows of 2 

9/30/2009 

1.39 

1539.4 


6 rows of 2 

9/30/2009 

1.14 

2032 


6 rows of 2 

9/30/2009 

1.46 


Bad Data 

6 rows of 2 

9/30/2009 

1.5 


Bad Data 

6 rows of 2 

9/30/2009 

1.46 

1036.7 


6 rows of 2 

9/30/2009 

1.41 

2032 


6 rows of 2 

9/30/2009 



Bad Data 

6 rows of 2 


10/2/2009 



Bad Data 

8 rows of 2 

10/2/2009 



Bad Data 

8 rows of 2 

10/2/2009 



Bad Data 

8 rows of 2 

10/2/2009 

1.13 


Bad Data 

8 rows of 2 

10/2/2009 

1.08 

2309.1 


8 rows of 2 

10/2/2009 

1.04 

2309 


8 rows of 2 

10/2/2009 



No Data 

8 rows of 2 

10/2/2009 

1.01 

2540 


8 rows of 2 

10/2/2009 

0.95 

2309 


8 rows of 2 

10/2/2009 



Bad Data 

8 rows of 2 

10/2/2009 

1.68 

2674 


8 rows of 2 

10/2/2009 

1.61 

4233.3 


8 rows of 2 

10/2/2009 

1.52 

1539.4 


8 rows of 2 

10/2/2009 

1.45 

8466.7 


8 rows of 2 

10/2/2009 



Bad Data 

8 rows of 2 

10/2/2009 

1.64 

2419.05 


5 rows of 2 and 3 inline on the center 

10/2/2009 



Bad Data 

5 rows of 2 and 3 inline on the center 
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Date 

C|J 

Wavespeed 

(m/s) 


Obstacle Configuration 

10/6/2009 

1.64 

2674 


5 rows of 2 and 3 inline on the center 

10/6/2009 

1.77 

2419 


5 rows of 2 and 3 inline on the center 

10/6/2009 



Bad Data 

5 rows of 2 and 3 inline on the center 

10/6/2009 

1.53 

2540 


5 rows of 2 and 3 inline on the center 

10/6/2009 

1.46 


Bad Data 

5 rows of 2 and 3 inline on the center 

10/6/2009 



Bad Data 

5 rows of 2 and 3 inline on the center 

10/6/2009 



Bad Data 

5 rows of 2 and 3 inline on the center 

10/6/2009 

1.37 

2540 


5 rows of 2 and 3 inline on the center 

10/6/2009 

1.27 

2540 


5 rows of 2 and 3 inline on the center 

10/6/2009 

1.17 

2419 


5 rows of 2 and 3 inline on the center 

10/6/2009 



Bad Data 

5 rows of 2 and 3 inline on the center 

10/6/2009 

1.04 

2309 


5 rows of 2 and 3 inline on the center 

10/6/2009 

0.94 

2117 


5 rows of 2 and 3 inline on the center 

10/6/2009 



Bad Data 

5 rows of 2 and 3 inline on the center 


10/6/2009 

1.72 

2309 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.71 

2309 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.62 

2419 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.51 

2419 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.43 

2674 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.34 

2540 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.25 

2674 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.16 

2419 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.07 

2419 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.06 

2309 


4 rows of 2 and 3 inline on the center 

10/6/2009 

1.01 

2309 


4 rows of 2 and 3 inline on the center 

10/6/2009 

0.9 

1881 


4 rows of 2 and 3 inline on the center 

10/6/2009 

0.8 



4 rows of 2 and 3 inline on the center 

10/6/2009 

1.71 

2209 


3 rows of 2 and 3 inline on the center 

10/6/2009 



Bad Data 

3 rows of 2 and 3 inline on the center 

10/6/2009 

1.6 

2117 


3 rows of 2 and 3 inline on the center 

10/6/2009 

1.53 

2419 


3 rows of 2 and 3 inline on the center 

10/6/2009 

1.51 


Bad Data 

3 rows of 2 and 3 inline on the center 

10/6/2009 

1.41 

2309 


3 rows of 2 and 3 inline on the center 

10/6/2009 



Bad Data 

3 rows of 2 and 3 inline on the center 

10/6/2009 

1.42 

2419 


3 rows of 2 and 3 inline on the center 

10/6/2009 

1.34 

2540 


3 rows of 2 and 3 inline on the center 

10/6/2009 

1.25 

2419 


3 rows of 2 and 3 inline on the center 

10/6/2009 

1.17 

2540 


3 rows of 2 and 3 inline on the center 

10/6/2009 

1.09 

2117 


3 rows of 2 and 3 inline on the center 

10/6/2009 

0.99 

2032 


3 rows of 2 and 3 inline on the center 

10/6/2009 



Bad Data 

3 rows of 2 and 3 inline on the center 
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Date 

C|J 

Wavespeed 

(m/s) 


Obstacle Configuration 

10/8/2009 

1.76 

2309 


6 Rows Inline Top/Bottom 

10/8/2009 

1.68 

2117 


6 Rows Inline Top/Bottom 

10/8/2009 

1.6 

2419 


6 Rows Inline Top/Bottom 

10/8/2009 

1.5 

2309 


6 Rows Inline Top/Bottom 

10/8/2009 

1.42 

2209 


6 Rows Inline Top/Bottom 

10/8/2009 

1.32 

2419 


6 Rows Inline Top/Bottom 

10/8/2009 



Bad 

Data 

6 Rows Inline Top/Bottom 

10/8/2009 

1.64 

2117 


5 Rows Inline Top/Bottom 

10/8/2009 

1.73 

2309 


5 Rows Inline Top/Bottom 

10/8/2009 

1.55 

2209 


5 Rows Inline Top/Bottom 

10/8/2009 

1.46 


Bad 

Data 

5 Rows Inline Top/Bottom 

10/8/2009 

1.37 

2309 


5 Rows Inline Top/Bottom 

10/8/2009 

1.3 

2117 


5 Rows Inline Top/Bottom 

10/8/2009 

1.23 

2309 


5 Rows Inline Top/Bottom 


** Bad Data Was Due to Faulty Kistler Pressure Transducers 
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